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We investigate Bianchi I cosmological model in the theory of a dilaton field coupled to gravity through a 
Gauss-Bonnet term. Two type of cosmological singularity are distinguished. The former is analogous to the 
Einstein gravity singularity, the latter (which does not appear in classical General Relativity) occurs when the 
main determinant of the system of field equations vanishes. An analogy between the latter cosmological singularity 
and the singularity inside a black hole with a dilatonic hair is discussed. Initial conditions, leading to these two 
types of cosmological singularity are found via numerical integration of the equation of motion. 
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General Relativity (GR) as leading theory of 
gravitational interaction at the classical level 
demonstrates a lot of interesting properties, for 
example, the existence of space time singularities 
|l]. The main "negative" feature of GR is that 
it is not renormalisable, so, it can not be quan- 
tized by the usual quantum field theory methods. 
The possible way to construct the quantum the- 
ory of gravity is to use string/M theory as the 
most promising candidate for the "theory of all 
physical interactions". After compactification to 
four dimensional space time one obtains some ef- 
fective gravitational theory which includes GR as 
the zeroth perturbation. This approach contains 
more wide set of solutions, and, hence, in addi- 
tion to the singularity classes introduced by GR 
can include new types of singular behavior. 

The simplest version of effective four dimen- 
sional string gravity action with the higher order 
curvature corrections up to the second order con- 
tains a coupling of the dilaton with gravity via 
the Gauss-Bonnet term as follows (we use the 
units where topl/n/Stt = 1): 
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where R is Ricci scalar, <j) is dilaton, A is dilatonic 
string coupling constant which is positive. Second 



order curvature correction Xe"'^'^ Sgb represents 
the product of the coupling functions e^^*^ and 
Gauss-Bonnet combination Sqe = RijkiR^-''''' — 
ARijR^^ + R^ . The solutions of this model have 
extensively been studied in the literature, both in 
a perturbative |^ and numerical |j,^ approaches. 
One important aspect is that when the second 
order curvature correction is taken into account a 
wide class of solutions contains new types of sin- 
gularities (earlier it was studied in d > 4 space 
time case, see, for instance, [|| and references 
therein). For example, during the investigation 
of the black hole solution (asymptotically flat, 
spherically symmetric and static) a new type of 
black hole inner singularity was found. It occurs 
at the finite (non zero) radius inside black hole 
and has the topology S^ x R^ (an infinite tube 
in time direction). This is curvature singularity, 
though it is week in Tipler's terminology m. If 
one works in curvature gauge 



ds^ = Adi^ 
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where A = A(r) and a = cr(r), Einstein- 
Lagrange equations can be written in a matrix 
form 



Oil A" -t- ai2(j' + Oizcf)" = bi, 
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where i = 1,2,3, a.y — a^ (A, A', ct, 0, (/)') and 
b,^b,{A,A',a,(t>A')- 



The describing singularity (we will call it as 
"determinant singularity") occurs when the sec- 
ond part (in quadratic brackets) of the main sys- 
tem determinant 



AA^ +BA + C 
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{A, B and C depend upon A', cr, 0, 0') vanishes. 
The asymptotic behavior of the metric and dila- 
tonic functions near this singularity r^ is H 
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where di, (T^, (\)i are numerical factors and 
r — Ts ^ 1. 

A numerical example is shown on Fig.|l|(a,b) 
taken from B. The solution was obtained with 
the help of the numerical method based on the in- 
tegration over an additional parameter along the 
solution trajectory described in Ref. |^. Some 
new mathematical aspects of this strategy were 
additionally studied in Ref Q . The solution con- 
sists of two branches which merge at determinant 
singularity Vs and solution could not be contin- 
ued further in radial coordinate. It is necessary to 
point out that all the metric functions (indepen- 
dently of their maximal derivative order) work as 
const/-yr — Ts in Riemannian tensor, therefore, 
cutvature square / = Rijutd'^^^ diverges. 

The main goal of this paper is to indicate that 
this kind of singularity is also significant in cos- 
mology. Previously, it was found in rather compli- 
cated models such as those elaborated to describe 
a dynamical compactification (we think that it 
is the first appearance of this kind of singularity 
in cosmology) [[lO| , in multidimensional Lovelock 
gravity |ll| and in Bianchi I |l2yi3t| and Bianchi 
IX |13| with moduli fields. In our paper we show 
that the determinant singularity occurs in the 
Bianchi I cosmology (the line element is 

ds^ = dt^ - a^(t)dxl - h^{t)dxl - c^{t)dxl, (5) 

a{t), b{t), c{t) are the scale factors) in the theory 
with the simplest stringy motivated second order 
curvature correction (1). 



Introducing three Hubble parameters p(i), q{t), 
r(t) according to following definition 

p(t) = -, qit) = l, r(t) = -, 
a c 

where the dot denotes d /dt, the variation of ac- 
tion (|l|) with respect to the metric components 
and the dilaton field gives the following equations 
of motion: 

pq + qr + rp + 2Afpqr - 0^/2 = 0, (6) 

(l-f 8r/)(g + 92) + (1 + 89/)(f + r2) 

+ (1 + %f)qr + (/.V2 = 0, (7) 

(1 + 8r/)(j9 + p2) + (1 + 8p/)(r + r^) 

+ (l + 8/)rp + ,/.V2 = 0, (8) 

(l + 8(?/)(p + p2) + (i + 8p/)(9 + g2) 

+ (1 + 8/)m + 072 = 0, (9) 

cf) ^ —{p + q + r)cj) + 8f'{pqr + pqr 

+ pqr + pqr{p + q + r)), (10) 

where /(i) — j^e~^'^, and the prime denotes dif- 
ferentiation with respect to (j). For mostly obvious 
presentation of the numerical results one have to 
introduce new variables h{t), a{t), (3{t) 

p = /i + a + V3/3, g = /i + a — v3/3, r = h — 2a, 

where /i is a generalization of a Hubble parameter 
of the isotropic case and describes the expansion 
rate of the universe, a and f3 are parameters of 
anisotropy. 

It should be emphasized that the presence of 
the Gauss-Bonnet term can in principle allow the 
violation of the weak energy condition B . This 
means that it is possible for the spatially flat uni- 
verse to make a transition from expansion to con- 
traction (see, for example, [g[) what can also be 
easily seen directly from the constraint equation 
(|^. Indeed, in Einstein gravity the geometrical 
part of constraint (|g) is 3h^ — 3a^ — 3/3^ which 
is nonpositive if /i = and, hence, this case is 
inconsistent with the constraint equation. Even 
if the Gauss-Bonnet term is taken into account in 
the isotropic case the correction 24,fpqr vanishes 
at ft. = 0, so the Einstein gravity result {h can not 
change it's sign in the spatially flat homogeneous 
universe) remains true. However, in Bianchi I 



case this correction contains a product of Hubble 
parameters which can be positive when the sum 
of the Hubble parameters {h = l/3(p + q + r)) 
becomes zero and the condition /i = no longer 
contradicts the constraint (g). 
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Figure 1. The behavior of the metric functions 
A (presents as d^A/dr^ in field equations) and 
a (presents as da / dr in field equations) against 
radial coordinate r in black hole case when the 
second order curvature corrections are taken into 
account 



The equations of motion (H)-(|o|) can be inte- 
grated numerically. They have 5 degrees of free- 
dom. Using constraint (^ one can reduce this 
number to 4. As it was done previously p| we 



calculate all dynamical variables checking the ac- 
curacy with the help of the constraint equation 
(g). In our calculation the left-hand side of the 
constraint equation did not exceed 10~^ unless 
we are very close to an ordinary singularity. In 
all our numerical investigations we put A = 1 for 
simplicity. 
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The determinant singularity occurs when Dmain 
vanishes. Our numerical investigations show that 
it can happen for a rather wide set of initial data 
(see below) . The asymptotic form of Hubble pa- 
rameters behavior near this singularity ts is (we 
checked it numerically and solved the equations 
on expansion coefficients /i^^, there are no any 
kinds of contradictions) 



hi = his + hii^/t -ts + hi2Wt - ts) +■■■, 
= (l)s + (l>i{V^^^^f + MV^^^sf +(,11) 

where hi = p^q,r and t — t^ <C 1. 

Comparing (11) with (4) we can see that these 
two expansions are absolutely analogous. The ex- 
pansion properties of functions having the same 
maximal derivative order in the equations of mo- 
tion (the first order for a and hi, the second or- 
der for A, and </>) are identical. It means that 
the properties of the determinant cosmological 
singularity are the same as the properties of the 
determinant Schwarzschild black hole singularity 
(keeping in mind, of course, that metric func- 
tions in cosmological case depend on time while in 
black hole they depend on the radial coordinate). 

A numerical example is plotted in Fig.0. 
Again, we can see two branches of solution and we 
can not continue these solutions further in time. 

In our investigation we calculated two- 
dimensional slices in (a, /?)-plane of initial condi- 
tion space with fixed values of h and (j). The ini- 





Figure 2. The behavior of the Hubble functions 
hi (presents as dhi/dt in field equations) and (f) 
(presents as d^cfi/dP in field equations) against 
coordinate t in Bianchi f case when the second 
order curvature corrections are taken into account 



tial value of (p has been taken from the constraint 
equation, which is quadratic with respect to (f>. 
The results below corresponds to the choise of 
the largest root of the quadratic equation. Each 
of the slices represent a 100 x 100 grid. For each 
node of the grid we calculated both future and 
past solutions over the time interval of 10^ in our 
units. There are three possible outcomes: ordi- 
nary singularity, determinant singularity and a 
nonsingular behavior. Such slices have been cal- 
culated for the values of (f> in the interval from 
6 — —15 to (h — +15 and for the values of h in 



the interval from h = 1 ■ 10~^ to ft. = 1. 

Each of the numerical solutions has been an- 
alyzed and classified. Finally, for each slice we 
obtained picture like one shown in Fig. y-|6|. Here 
the black color corresponds to non-singular solu- 
tions. The dark-gray color corresponds to ordi- 
nary singular solutions and the light-gray color 
corresponds to solutions with determinant singu- 
larity. The white color indicates prohibited re- 
gions in initial condition space where none of the 
solutions can start. 

For large values of </> the allowed region in 
(a, /3)-plane has a compact subregion of approx- 
imately circular shape with the center in point 
(0,0) (see Fig.H). This can be easily seen from 
constraint equation (o) if the term proportional 
to e"^"^ is neglected. Numerical calculations show 
that in this case all future solutions are non- 
singular while all past solutions are singular with 
Dmain ^ 0. There is only one exception from 
this rule: the past solution for a = /3 = (the 
isotropic case) has an ordinary singularity. 

With decreasing of (f) this picture changes dra- 
matically. The shape of allowed region undergoes 
substantial transformation (see Fig.g): the for- 
mer compact subregion becomes connected whith 
the outer non-compact one. The more important 
feature is that some trajectories meet a singu- 
larity in future, either due to a recoUapse (see 
above) or vanishing of the main determinant. The 
fraction of such solutions increases with initial (j) 
decreasing. For some value of initial the solu- 
tions, nonsingular in future, disappear completely 
(apart from the exceptional isotropic case). This 
critical value depends on generalized Hubble pa- 
rameter h as it is plotted in Fig. 0. The transi- 
tion from the regime with only nonsingular fu- 
ture solutions to one with only singular future 
solutions is rather sharp according to our results 
which can be easily seen from Fig. 0. The uni- 
verse with (j) < (f)Q (for a given h) and an arbitrary 
small nonzero initial anisotropy can not leave a 
high-curvature regime and falls into a singular- 
ity, so such initial conditions are not suitable for 
describing the evolution of the our Universe. 

The other important dynamical feature that 
was found via numerical integrations is related 
with past asymptotic. Though both described 
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Figure 3. The minimal initial values of cf) which 
allow solutions non-singular in future (solid line) 
and the maximal initial values of </> which allow 
singular solutions in future (dot-dashed line). 



types of singularity are valid for various initial 
data sets, all the initial data which lead to non- 
singular future behavior have determinant singu- 
larity in the past (again, we ignore the exceptional 
isotropic solution). In particular, for large initial 
(j) the only possible past singularity is the deter- 
minant singularity (Fig. ^ . 

To obtain this result we investigated a rather 
wide set of initial conditions. However, the full set 
of allowed initial conditions is noncompact and, 
hence, we have no proof of this statement in the 
strict sense. Hence, we can put forward a con- 
jecture: the past singularity for nonsingular in 
future Bianchi I cosmological model with the sec- 
ond order curvature corrections is of determinant 
type. This proposal is analogous to the result ob- 
tained in p] for Schwarzschild black hole, where 
the asymptotically Minkowski initial data at in- 
finity lead to a determinant singularity at the fi- 
nite values of the radial coordinate. 
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Figure 4. A slices of initial condition space in (a,/3)-plane for h = 0.16 and = 0. Figure (a) corresponds 
to the past solutions and figure (b) corresponds to the future solutions. The black color corresponds to 
non-singular solutions, the dark-gray color corresponds to ordinary singular solutions and the light-gray 
color corresponds to solutions with determinant singularity. The white color indicates prohibited regions 
in the initial condition space. 
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Figure 5. The similar slices for h = 0.16 and (p — —1. 
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Figure 6. The similar slices for h = 0.16 and = — 2. 



